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Aging dynamics in thin films of poly(methyl methacrylate) (PMMA) have been investigated 
through dielectric measurements for different types of aging processes. The dielectric constant was 
found to decrease with increasing aging time at an aging temperature in many cases. An increase 
in the dielectric constant was also observed in the long time region (>1 fh) near the glass transition 
temperature for thin films with thickness less than 26nm. In the constant rate mode including a 
temporary stop at a temperature T a , the memory of the aging at T a was found to be kept and then 
to be recalled during the subsequent heating process. In the negative temperature cycling process, a 
strong rejuvenation effect has been observed after the temperature shift from the initial temperature 
T\ to the second temperature ?2(= Ti + AT) when AT w — 20K. Furthermore, a full memory effect 
has also been observed for the temperature shift from Ti to Ti. This suggests that the aging at Ti is 
totally independent of that at T2 for AT ~ — 20K. As |AT| decreases, the independence of the aging 
between the two temperatures was found to be weaken, i.e., the effective time, which is a measure 
of the contribution of the aging at Ti to that at T2, is a decreasing function of |AT| in the negative 
region of AT. As the film thickness decreases from 5I4nm to 26nm, the |AT| dependence of the 
effective time was found to become much stronger. The contribution of the aging at T2 to that at 
Ti disappears more rapidly with increasing |AT| in thin film geometry than in the bulk state. 

PACS numbers: 71.55.Jv; 81.05.Lg; 77.22.Ch 



I. INTRODUCTION 

In glassy materials it is well-known that a very slow 
relaxation towards an equilibrium state is observed be- 
low the glass transition temperature T g 0, . This slow 
relaxation is called aging and is regarded as an important 
common property characteristic of disordered materials 
such as spin glasses 0,0, IE Ht > orientational glasses , 
supercooled liquids Tp, relaxor ferroelectrics Q , and 
polymer glasses 0, lllf . The investigations on such 
disordered systems revealed that memory and rejuvena- 
tion effects are observed during the aging process. The 
characteristic behavior of the aging dynamics below T g 
is closely related to the nature of the glass transition. 
Therefore, it is expected that the elucidation of the aging 
phenomena can lead to full understanding of the mecha- 
nism of the glass transition in disordered materials, which 
is still a controversial topics ^2 • 

In polymeric systems, Kovacs et al. investigated the 
memory effects observed in the aging phenomena |13|. 
In their experiments, a volume relaxation was measured 
on a glassy state located on the line extrapolated from 
the equilibrium line corresponding to a liquid state. Be- 
cause the glassy state is on the equilibrium line, it is 
expected that no further volume change occurs on this 
state. However, their results suggest that the volume re- 
laxation strongly depends on the thermal history which 
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the polymer had experienced before arriving at this equi- 
librium state. In other words, the thermal history can be 
memorized within the glassy state of polymers. Recent 
experiments by Miyamoto et al. revealed that polymeric 
materials can keep not only thermal history but also me- 
chanical history in their glassy states and the histories 
can be recalled upon the glass transition, through mea- 
surements of tension in rubbers |14| . 

In spin glasses, interesting studies on the aging phe- 
nomena have been done in recent years 0. El E3 In the 
experiments on the CdCri. 7^0.384 insulating spin glass, 
a relaxation of the imaginary part of the magnetic ac- 
susceptibility \" was observed in a temperature cycling. 
First, the sample is cooled down to a temperature Ti be- 
low T g and is kept at this temperature for a period of 
Ti . Then the temperature is changed from Tj to T\ + AT 
(AT < 0) and kept for T2, and after that the temperature 
is returned to T\ . It was observed that upon cooling from 
Ti to T\ + AT, the system rejuvenates and returns to a 
zero-age structure even after a long stay at the higher 
temperature T± . The interesting point is that the system 
has a perfect memory of the past thermal history, i.e., 
when heated back to T\ , the susceptibility x" agrees well 
with the value that \" bad reached just at the end of the 
first stage of the aging at T\. As for polymeric systems, 
a similar measurement on the aging dynamics has been 
done for the dielectric constant in poly(methyl methacry- 
late) (PMMA) by Bellon et al. to investigate the memory 
effects of the aging dynamics. Although a memory effect 
similar to that in spin glas ses has been observed also in 
PMMA in their studies [lfj, |ll| , more detailed investiga- 
tions would be required to extract an universal behavior 
of the aging dynamics in polymer glasses and to compare 
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the aging dynamics between polymeric systems and spin 
glasses. 

It is very natural to believe that the aging dynamics 
in the glassy state mentioned above are related to the 
nature of the glass transition. Recent investigations on 
the glass transition in thin film geometry show that there 
is a strong dependence of the glass transition tempera- 
ture and the dynamics on film thickness and the system 
size Ha, El* Eil • Therefore, it is expected that the aging 
dynamics depend on the film thickness |20j| , and through 
such investigations it would be possible to approach the 
basic mechanism of the glass transition. From the tech- 
nical point of view, thin film geometry is ideal for di- 
electric measurements because the smaller thickness of 
the films enhances the resolution of capacitance mea- 
surements. For this reason, precise measurements can 
be done for thin film geometry. 

In this study, we measured time evolution of the com- 
plex electric capacitance of PMMA films during the ag- 
ing process using dielectric relaxation spectroscopy. The 
purpose of this study is to investigate the relaxation dy- 
namics of non-equilibrium states including the memory 
and rejuvenation effects in polymeric systems and to elu- 
cidate the thickness dependence of the aging dynamics. 

This paper consists of six sections. After giving exper- 
imental details in Sec. II, the aging dynamics at a given 
temperature are shown in Sec. III. In Sec. IV, experimen- 
tal results on the memory and rejuvenation effects are 
presented for two different thermal histories. After dis- 
cussing the experimental results in Sec.V, a summary of 
this paper is given in Sec. VI. 



trolled through a heater wound around the cell using a 
temperature controller. The temperature of the poly- 
mer sample is monitored as the temperature measured 
through a thermocouple put onto the back side of the 
glass substrate. Heating above 400 K was done several 
times before the measurements so that the contact be- 
tween the polymer films and the substrate is stabilized. 
The temperature of the sample is stabilized within 2 min- 
utes after a temperature jump. The stability of the tem- 
perature control is ±0.IK. 

Dielectric measurements were done by an LCR meter 
(HP4284A). The frequency range of the applied electric 
field was from 20Hz to I MHz. The voltage level was con- 
trolled so that the strength of the applied electric field 
remains almost constant independent of film thickness. 
The typical electric field was 2xl0 6 V/m. In our mea- 
surements, the complex electric capacitance of the sample 
condenser C*(T) was measured and then converted into 
the dynamic (complex) dielectric constant e*(T) by di- 
viding the C*(T) by the geometrical capacitance Co(T ) 
at a standard temperature To. The value of C* is given 
by C* — £*£()■§ and C'o — eo4, where eo is the permit- 
tivity in vacuum, S is the area of the electrode and d is 
the film thickness. For evaluation of e* and Co, we use 
the thickness d which is determined at To =273 K and 
S=8x 10~ 6 mm 2 . Because e* and C* are complex num- 
bers, we define the real and imaginary parts as follows: 
e*=e' - it" and C*=C" - iC" . 



III. AGING DYNAMICS 



II. EXPERIMENT 

Polymer samples used in this study are atactic PMMA 
purchased from Scientific Polymer Products, Inc. The 
weight-averaged molecular weight M w =3.56xI0 5 and 
M w /M n =1.07, where M n is the number-averaged molec- 
ular weight. The glass transition temperature T g de- 
termined by differential scanning calorimetry (DSC) is 
about 380 K for the bulk samples of PMMA used in this 
study. The value of T g in thin films with <i=26nm is 
lower by 7K than T g in the bulk samples. Thin films 
are prepared on aluminum vacuum-deposited glass sub- 
strate using a spin coat method from a toluene solution 
of PMMA. The thickness of the films is controlled to be 
20nm to 514nm by changing the concentration of the so- 
lution of PMMA. The thickness is evaluated from the 
value of the electric capacitance at 273 K of as prepared 
films before measurements in the same way as in previous 
papers ^10,^3- After annealing at 343K, aluminum 
is vacuum deposited again to serve as upper electrodes. 
A copper wire is attached to the end of the electrode 
with conducting paste so that the polymer samples can 
easily be connected with the measurement system. The 
sample prepared in the above way is mounted inside a 
sample cell. The temperature of the sample cell is con- 



In this section, we show the relaxation behavior during 
the isothermal aging process at a given temperature T a . 
In order to remove the thermal history of the sample, the 
system is heated up to 403K (> T g ) and is kept at 403K 
for a couple of hours. After this procedure the system is 
cooled down to T a at the rate of 0.5 K/min and the mea- 
surement of the complex electric capacitance is started. 
Figure 1(a) displays the time evolution of the imaginary 
component of the dynamic dielectric constant e" at vari- 
ous aging temperatures for the thin films with thickness 
of 20nm. In Fig. I (a), it is found that e" is a monotonical 
decreasing function of the aging time t w . The vertical 
axis in this figure is the value of e" normalized with re- 
spect to the initial value e"(0). The origin of the time 
axis is defined as the time at which the temperature of 
the system reached T a . The decrease in the dielectric 
constant corresponds to the smaller response against an 
applied external field, and hence this indicates that the 
system goes down to a deeper valley of the free energy 
via thermal activations within the energy landscape pic- 
ture |2ll |22| . In other words, the system approaches an 
equilibrium state as the time t w elapses at a temperature 
T a below T g . The aging phenomena can be interpreted 
as this approach to the equilibrium state. 

In an equilibrium state, the dielectric constant can be 
uniquely specified as functions of temperature T, pres- 
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FIG. 1: Aging time dependence of the dynamic dielectric 
constant normalized with the value at t w =0 in PMMA thin 
films with d=20nm for various aging temperatures: T a =375.4 
K, 365.6 K, 355.8 K, 346.0 K, and 336.1 K: (a) the imaginary 
component and (b) the real component. The frequency of the 
applied electric field / is 20Hz. The horizontal axis is the 
natural logarithm of t w . The solid curves are obtained by 
fitting the data points with Eq.(5). The best-fit parameters 
for the imaginary parts are listed in Table I. 



TABLE I: Fitting parameters obtained by fitting the ob- 
served values of e" with Eq.(5) for thin films of PMMA with 
d = 20nm. The frequency of applied electric field is 20Hz. 
The exponent n is fixed to be 0.309, which is obtained by 
fitting the data at 355.8 K 



Ta (K) 


Ae' a ' 9 A"(0)(xl0- 2 ) 


to (min) 


£ "(0) 




375.4 


9.286 ±0.007 


21.8±0.1 


0.1147 


0.01065 


365.6 


5.486 ± 0.007 


18.6±0.2 


0.1302 


0.00714 


355.8 


3.081 ±0.041 


26.4±0.9 


0.1601 


0.00493 


346.0 


1.989 ± 0.006 


52.3±0.7 


0.2016 


0.00401 


336.1 


1.357 ±0.007 


35.3±0.9 


0.2495 


0.00339 



sure p and the frequency of the applied electric field / 
for the polymer films with a fixed thickness. For simplic- 
ity, we ignore the effect of the pressure in the discussions 
below. If the system is in a non-equilibrium state, it is 
not easy to specify the dielectric constant, and the aging 
phenomena must be taken into account. The dielectric 
constant depends on the thermal history that the sample 
has experienced. In general, if the system is aged at a 



temperature T Qi i for a period of t w> i, at T aj2 for t Wy2 , ■ ■ ■ , 
and at T a>n for t w _ n , the dielectric loss e" can be specified 
by the following description : 



= e"(T,u;{t wd (T ad )}) ) 



(1) 



where u = 2nf and {t wJ (T a ^)} = {t wA (T aA ),t w ^(T^ 2 ), 
■ ■ ■ t w>n (T a>n )}. If the system is aged for t w at just one 
temperature T a , the dielectric loss at the temperature T 
is described by 



e"(T,w;t w (T a )). 



(2) 



Firstly, the aging dynamics under an isothermal condi- 
tion at T=T a are considered. In principle, the dielectric 
loss e"(T Q , w; oo) is the value of e" in the equilibrium state 
at T a and the rest of e" can be defined as the aging part 
e''(T a ,u>;t w ) in the following way: 



'(T a ,u,t w ) = e"(T a ,w;oo) ± e''{T a ,u);t w ). 



(3) 



The symbols T a and u>, here, are dropped for simplicity 
and it is assumed that the aging part of e" is given by 
the following equation: 



£ oo(*i») ~~ 



"9 



(i ± t w /t ) r 



(4) 



where Ae" g is the relaxation strength towards the equi- 
librium value, to is the characteristic time of the aging 
dynamics, and n is an exponent. If t w is much larger than 
to, e" g ~ t~ n . In this case, the dielectric loss e" decays 
according to the power-law decay in the long time region. 
The similar expressions for e' are also used in this paper. 

The solid curves given in Fig. 1(a) were obtained by 
a non-linear least square fit of the observed data to the 
following equation derived from Eqs.(3) and (4): 



e"(0) 



A4' g 
e"(0) 



(5) 



Table I displays the fitting parameters obtained by Eq. JSJ) 
for e" in thin films of PMMA with d=20nm at various 
aging temperatures, where the value of n is fixed to be 
the value obtained for 355.8 K. The relaxation strength 
Ae" g was found to decrease with decreasing aging tem- 
perature T a . Furthermore, the characteristic time to of 
the relaxation of e" in the aging process was found to 
increase with decreasing T a , in other words, the aging 
dynamics become slower with decreasing T a . Because 
the observed data in Fig. 1(a) can well be reproduced us- 
ing Eq.(4), it is implied that the relaxation of e" at an 
isothermal condition below T g obeys the power-law one 
for the long-time region in the thin films with d = 20nm. 
According to the report in Ref. 10], the aging behav- 
ior depends on the frequency of the applied electric field. 
As the frequency decreases, the relaxation strength due 
to aging increases for the frequency range from 0.2Hz to 
20Hz. In this study, we have measured the dynamic di- 
electric constant over 4 decades and could confirm the 
results in Ref. [ToT |. 
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there is no tendency to saturate to a constant value. In 
this case, for t w > 20min, the relaxation of e" due to ag- 
ing obeys a logarithmic law. This logarithmic decay law 
was observed also in the measurements done by Bcllon 
et al. Although it is impossible to judge whether the 
logarithmic decay is valid also outside the time window, 
there may be a qualitative difference between the decay 
law for the thin films and that for the thicker films. 

In Sec. Ill, we have observed that there is the de- 
crease in e' and e" with increasing t w at various values of 
T a (— 336. 1K~ 375. 4K), except under the condition that 
t w > llh, d <26nm, and T a =374K, where the increase 
in e' with t w was observed. In the discussion section, the 
physical origin of the time variation of e' and e" will be 
discussed. 



IV. MEMORY AND REJUVENATION EFFECTS 



FIG. 2: Aging time dependence of the dielectric loss as a 
function of film thickness for PMMA thin films: d=20nm, 
26nm, 78nm and 514nm, T a =375K and /=20Hz. A deviation 
from the straight line implies that the logarithmic decay does 
not hold in the long time region for d=20nm and 26nm. 



As mentioned later, under the present experimental 
conditions that T a = 336. 1K~ 375.4K and / = 20Hz are 
located between the a-process and the /3-process in the 
dispersion map. Because Ae" g decreases with decreasing 
T a , the relaxation of e" due to the aging may be associ- 
ated with the a-process. 

Figure 1 (b) displays the t w dependence of the real part 
of the complex dielectric constant e' normalized with re- 
spect to the initial value at various aging temperatures 
T a . In this figure, it is found that e' decreases with in- 
creasing t w in a similar way to e" for T a 's less than 365.6 
K. However, there is a surprising difference between e 
and e" for T a =375 K: e' decreases with t w up to t w w 
llh and then begins to increase with t w . This increase in 
e' in the long time region was not observed for lower val- 
ues of T a . A similar increase in e' was observed also for 
d =26 nm. On the contrary, for thicker films (d =78, 514 
nm), no such increase in e' was observed in the present 
measurements. 

Wc here move to the results for the thin films with var- 
ious film thickness to see the thickness dependence of the 
aging dynamics. Figure 2 displays the time variation of 
the dielectric loss e" at T a =375 K for the thin films with 
thickness between 20nm and 514nm. Because the hori- 
zontal axis in Fig. 2 is on a logarithmic scale, any straight 
line indicates the existence of a logarithmic decay of the 
dielectric loss. For d=20nm and 26nm, we can see that 
there is a slight deviation from the straight line in the 
short and long time regions. This is consistent with the 
result that the time variation in e" obeys Eq.(5) for d=20 
nm, as shown in Fig. 1(a). For d =78nm and 514nm, on 
the other hand, e" decreases linearly with logt w for the 
region of t w > 20 min. Within the present time window 



In Sec. IV we show the memory and rejuvenation ef- 
fects observed during the aging process and discuss the 
detailed properties including the thickness dependence. 
For this purpose, we adopt two modes of thermal his- 
tory: 1) constant rate mode (CR mode) - After the ther- 
mal equilibration at 403 K, the sample is cooled at a 
rate of 0.5 K/min down to a temperature T a , and then 
is isothermally annealed (aged) at T a for 10 hours. The 
sample is then cooled down to room temperature at 0.5 
K/min. Finally, the sample is heated again from room 
temperature to 403 K at 0.5 K/min. 2) temperature cy- 
cling mode (TC mode) - The sample is cooled at a rate 
of 0.5 K/min from 403 K down to a temperature T\ and 
is aged at T\ for a period of tj.. Then, the temperature 
is rapidly changed from T x to T 2 (T 2 = T x + AT) and 
is kept at T 2 for t 2 . Finally, the temperature is changed 
back from T 2 to T\ and is kept at T\ for T3. In this 
section, we show the data for t\ — t 2 — 5h. Both the 
modes have originally been developed in order to investi- 
gate the relaxation behavior in the aging phenomena in 
spin glasses 0, 0, Q . 

A. Constant rate mode 

As it is well-known, there arc two typical molecular 
motions in PMMA, which are usually called the a-process 
and the /3-process |23( • The a-process is strongly related 
to the glass transition, while the /3-process is due to the 
local motion of polymer chains. In the case of PMMA, 
the dielectric relaxation strength of the /3-process is larger 
than that of the a-process due to the existence of polar 
bulky side groups. Both the dielectric constant and loss 
have a strong temperature dependence due to the a- and 
/3-processes. For example, e" exhibits two distinct loss 
peaks in a temperature domain at a given frequency : 
one is at about 310 K due to the /3-process and the other 
is at about 400 K due to the a-process for / = 20Hz and 
d = 20nm. (See the inset of Fig. 6.) 
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FIG. 3: Temperature dependence of the real part of the com- 
plex electric capacitance e'Co observed by the CR mode with 
T a = 374. 2K for PMMA films with d=514nm. The frequency 
of the applied electric field / is 20 Hz. The symbols o and A 
stand for the cooling and heating processes of the CR mode, 
respectively. The solid curve (— ) and dotted curve (• • • ) are 
for the reference data observed during the cooling and heating 
processes without any temporary stop at T a . 



FIG. 4: Temperature dependence of the difference between e 
and e' le{ observed by CR mode with three different aging tem- 
peratures T a =374.2 K, 364.5 K and 354.5 K for PMMA films 
with d=514nm. The frequency of the applied electric field 
/ is 20Hz. The heating and cooling rates were 0.5 K/min 
and the aging times at T a was 10 hours. Open and full sym- 
bols display the data observed during the cooling and heating 
processes, respectively. 



For this CR mode, we measured both the dielectric 
constant e' and the dielectric loss e". If the range of t w 
is restricted to t w < 10 hours, there is no qualitative dif- 
ference in the aging dynamics between e' and e". In the 
case of <i=514nm, the errors attached to the value of e" 
are larger than those of e' because of the resolution of the 
LCR meter. For this reason, we show here the results of 
e', which can be described as e'(T, u>; t w (T a )) according 
to our notation. In addition to the e', we prepare a refer- 
ence curve, e'(T, ui;0), both for the heating and cooling 
processes at the same rate (0.5K/min) without any tem- 
porary stop during cooling process. We refer to the curve 
as4 f (=e'(T>;0)). 

Figure 3 shows an example of the temperature depen- 
dence of the real part of the complex electric capacitance 
e'C observed by the CR mode with T Q =374.2 K for the 
films with ri=514nm. The reference curves eJ. ef Co for the 
cooling process and the subsequent heating process are 
given by the solid and dotted curves, respectively. The 
curve of e(. ef during the cooling process is slightly dif- 
ferent from that during the heating process. The data 
points denoted by circles stand for e'Co observed during 
the cooling process, while those denoted by triangles do 
for e'Co observed during the heating process. 

Figure 4 displays the temperature change in e' — e' f , 
i.e., the difference between e' and e' le{ obtained by the CR 
mode with various aging temperatures T a for d = 514nm. 
Here, we use e' lef of the cooling (heating) process for e' 
of the cooling (heating) process. In Fig. 4, it is found 
that, as the sample is cooled down from 403 K to T a 
(=374. 2K), e' is almost equal to e' re{ and the difference 
e' — e ref remains zero. The difference then begins to devi- 
ate from to a negative value and the deviation becomes 
larger monotonically with increasing aging time (at least 
up to t w ~ 10 hours) during the subsequent isothermal 



aging at T a . This indicates that e' decreases with ag- 
ing time at T a . As the temperature decreases from T a 
to room temperature after the isothermal aging, the dif- 
ference e' — e' Ic{ decreases with decreasing temperature 
and becomes almost zero at room temperature. Just af- 
ter the temperature reaches 293 K, the heating process 
at the rate of 0.5K/min starts. As the temperature in- 
creases, the difference e' — e' rcf begins to deviate from to 
a negative value almost along the same way as observed 
during the preceding cooling process. The difference be- 
tween e' and e' Te{ exhibits a maximum at around T a +10K 
and then decreases to with increasing temperature. 

This behavior observed by the CR mode can be inter- 
preted as follows: the thermal history that the sample is 
aged at T a for 10 hours is memorized at T a during the 
cooling process. As the temperature decreases from T a to 
room temperature, the sample begins to rejuvenate and 
is back to almost a standard age at room temperature. 
During the subsequent heating process, the sample be- 
comes older almost according to the curve along which 
the sample experienced rejuvenation during the preced- 
ing cooling process after a temporary stop at T a . This 
result implies that not only the aging process at T a near 
T g but also the subsequent cooling process can be memo- 
rized and the whole thermal history can be read out dur- 
ing the heating process. A similar temperature change 
in e' — e' f has already been observed in PMMA films of 
thickness of 0.3 mm [lOj. 

To investigate the thickness dependence of this aging 
behavior, we show the results for various film thicknesses, 
d =26, 78, and 514nm. In Fig. 5, it is found that the 
strength of aging for 10 hours at 374 K decreases with 
decreasing film thickness. Furthermore, during the sub- 
sequent cooling process, the difference between e' and ej. cf 
approaches zero more quickly in the thin films than in the 
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FIG. 6: Temperature dependence of the dielectric loss ob- 
served by the TC mode with Ti= 375. 3K and T 2 = 355.8K 
(AT = -19.5K) for PMMA thin films with d=20nm and the 
frequency of the applied electric field /=20Hz. The dotted 
curve displays the temperature dependence of e" obtained 
during the cooling process with the rate of 0.5K/min without 
any isothermal aging. The overall behavior of e" for temper- 
ature range from 310K to 403K is shown in the inset. 



FIG. 5: Difference between e' and e' ref as functions of temper- 
ature observed by the CR mode for PMMA thin films with 
d=26nm, 78nm, and 514nm. The upper figure displays the 
results for the cooling process and the lower one does those 
for the subsequent heating process. Dielectric measurements 
were done for /=20Hz. The isothermal aging was done at 
T a = 374K during the cooling process. 

thicker films. In other words, the thin films are found to 
age more slowly, but to rejuvenate more quickly. 

B. Temperature cycling with negative AT 

In this subsection, we show the results obtained by the 
TC mode. In this mode, the sample is aged at two dif- 
ferent temperatures T± and T 2 , and hence e" is described 
as the quantity e"(T,u>] t w \, t W 2), where t w \ and t W 2 are 
the aging times at temperatures T\ and T 2 , respectively. 
Here, the total aging time t w is given by t w — t w \ + t w2 - 
Using this expression, the dielectric loss e" observed by 
the TC mode can be described in the following way: 

r e"(Ti,w; *„,,()) : < t w < n 
e" = < e"(T 2 , u; n, t w - n) : n < t w < n + r 2 (6) 
[ e"{Ti,w;t w - T2,t 2 ) : t w >T 1 +r 2 . 

For example, for a value of t w (ri < t w < n + r 2 ), the 
notation e"(T2, cu; n, t w — n) corresponds to the imagi- 
nary part of the dielectric constant observed by frequency 
/ = cu/2it at T 2 after the aging at T\ for n and at T 2 for 
t w —T\. In this case, the reference curve e" ef can also be 
defined as e^ f = e"(T,co; 0, 0). 



Figure 6 displays a typical behavior of e" and e" ef of the 
TC mode for PMMA thin films with d = 20nm and / = 
20Hz. In the inset, the temperature dependence of e" ef at 
20Hz is shown for the temperature range from 290K to 
403K during the cooling process (See the dotted curve). 
The two distinct loss peaks, the a-process and the j3- 
process, are observed. As the sample is cooled from 403K 
to T\(=375.3K) at the rate of 0.5K/min according to the 
procedure of the TC mode, the dielectric loss e" moves 
along the reference curve e" ef and then deviates from the 
reference curve to a lower value upon the beginning of 
the aging at T\ (See the solid curve). The main part of 
Fig. 6 shows this behavior in the temperature range from 
350 K to 380 K. As the temperature decreases from T\ to 
T 2 (=355.8K) by AT=-19.5 K after the aging at T x for a 
period of t\ , the dielectric loss e" immediately returns to 
a value on the reference curve at T 2 . As the temperature 
is kept at T 2 , e" deviates from the value on the reference 
curve. While the temperature is changed from T 2 to T\ 
and is kept at T\ after the aging at T 2 for a period of t 2 , 
e" goes back to a value at T\, which deviates from the 
reference curve, and then e" decreases with time during 
the aging process at T\. 

To elucidate the aging dynamics at a temperature in 
the TC mode more in detail, it is useful to investigate 
the difference between e" and e" ef as we have already 
done in the CR mode. In this case, the reference val- 
ues e" ef (Tj,w) at Tj (j = 1,2) are subtracted from the 
observed dielectric losses e" at Tj, respectively. Here, we 
define e" - e? ef = e"(T, w; t wl ,t w2 ) - e" ef (T, u>; 0, 0), where 
t w =t w \ + t W 2- For the TC mode, 
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FIG. 7: (a) Difference between e" and e" et observed by the TC mode with Ti=375.3K and T 2 =355.8K (AT=-19.5K) for 
PMMA thin films with d=20nm. (b) The difference e" — e" ei obtained by shifting the date points in the third stage in the 
negative direction of the time axis by t 2 after removing the data points in the second stage. Aging times at the first and second 
stages are n = T2 = 5 hours. The horizontal axis of Fig. 7(b) is the total aging time at Ti. The solid ( — ) and dotted curves 
(■ • • ) are standard relaxation ones obtained by isothermal aging at Ti =375.3 K and T2=355.8 K, respectively. It should be 
noted that the time origin of the dotted curve for Tb is shifted from t w =0 to t w =Ti. 



4 f ee e"(T,uj;t wU t w2 )-e'^(T,ij;0,0) 

e"{T 1) cj;t W) 0) - 4' ef (T\, w; 0, 0) 
= { €"{T 2 ,u;T U t w -T 1 )-e'; ei {T 2 ,^,Q) 
e"(Ti ,u;t w - t 2 ,t 2 ) - 4 f (Ti , uj; 0, 0) 



o < t w < n 

Tl < t w < Tl 
tw > Tl+T 2 



- T 2 



(7) 



Figure 7 displays the time dependence of e" — e" of for 
the thin films with d=20nm, Ti=375.3K, T 2 =355.8K, 
and AT=-19.5K. The origin of the time (t w =0) is de- 
fined as the time when the temperature reaches Ti after 
the cooling from 403 K (above T g ). The value of e" — e" cf 
is found to decrease with increasing t w and to change 
immediately to when the temperature is lowered from 
Ti to T 2 at t u , — ti. This indicates that the decrease 
in temperature from Ti to T 2 causes the polymer sam- 
ple to rejuvenate. As the aging time elapses at T 2 af- 
ter the temperature shift, e" — e" ef is found to exhibit 
a relaxation just like a new relaxation process starts at 
t w = ti. The dotted curve in Fig. 7(a) which starts at 
t w = Ti is the relaxation curve of e" — e" cf observed at T 2 
after the temperature change directly from 403 K to T 2 . 
The origin of the time axis is shifted from to t w = T\ 
for the dotted curve in Fig. 7(a). The observed data for 
Ti < t w < Ti + t 2 are located exactly upon this dotted 
curve. At t w = t\ + r 2 , the system is heated up to T\, 
and then the temperature of the system is kept at T\. 
On this stage, e" — e" cf is found to go back to the value 
which e" — e" of had reached at t w = t\ and to begin to 
decrease as if there were no temperature change during 
the aging process at Ti. 

In Fig. 7(b) is shown the time evolution of e" — e" ef after 
removing the data between ti and n + t 2 and shifting 
the data for t w > T\ + t 2 in the negative direction along 
the time axis by t 2 . We can consider the horizontal axis 



as the total time which the system spent at Ti. In this 
figure, we can see that e" — e" of decreases monotonically 
with the aging time without any discontinuous change 
except for a short region just after t w = T\. The curve 
obtained in the above way agrees very well with that 
obtained by keeping the system at T\ without any tem- 
perature change, e"(Ti, w\ t w , 0) — e" e{ (T^ w; 0, 0) (We call 
this curve the standard relaxation curve.). This implies 
that polymer glasses can remember the state at t w — ti 
of the relaxation towards the equilibrium state and re- 
call the memory at t w — t\ + t 2 . We call this behavior a 
complete memory effect. This was a result obtained for 
AT = -19.5K. 

We here show the aging behavior for various values 
of AT for PMMA thin films with d = 26nm. Fig- 
ures 8(a)-(d) show e" — e" of as a function of the aging 
time t w for various values of AT— —19.6 K to —4.8 
K and Ti =374. 2K. Figures 8(e)-(h) are obtained from 
Figs.8(a)-(d) according to the following procedure: the 
data points for t w = n ~ T\ + t 2 are removed and the 
ones for t w > ti + t 2 are shifted by an amount of At in 
the negative direction along the t w axis so that the data 
can be overlapped with the standard relaxation curve for 
the aging at T± in the most reasonable way. Here, an 
effective time T e g is defined as the time t 2 — At. The 
effective time can be considered as a measure of the con- 
tributions of the aging at T 2 to that at Ti. If T e ff is zero, 
the aging at T 2 has no contributions to that at T±, while 
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FIG. 8: Difference between e" and e" ef observed by the TC mode for PMMA thin films with d=26nm and /=20Hz. The aging 
temperature of the first stage Ti is fixed to be 374. 2K and T 2 (= Ti + AT) is changed: AT=-4.8 K (d,h), -9.9 K (c,g), -15.0 
K (b,f), and —19.6 K (a,e). Aging times are ti = ti = 5 hours. In Figs. 8(e)-(h), the data points of the second stages are 
removed and those of the third stage are shifted by an amount of time in the negative direction of the time axis so that the 
data points of the third stages can be overlapped with the standard relaxation curve most reasonably. The effective time r e g 
is evaluated as the difference between the time at the end of the first stage and the time at the beginning of the initial time of 
the third stage after the above data fitting process. The difference in length between the arrow of T2 in (a) and the arrow of 
r e ff in (h) shows how the aging at T 2 affects the aging at Ti . 



if T c ff is equal to T2, the aging at T 2 has full contributions 
to that at T\. In Figs.8(a)-(d), it is found that the value 
of e" - e" ef is reinitialized at t w = n for AT=-19.6K~ 
— 9.9K. However, the reinitialization at t w = ti is not 
complete for AT=-4.8 K. 

As shown in Figs.8(e)-(h), r e ff=0 for 
AT=-19.6K 9.9K, and r off w 200min for AT=-4.8 



K. This indicates that for |AT| >9.9K there is a 
complete memory effect and the aging at T 2 gives no 
contributions to that at Ti, in other words, the aging 
at Ti is independent of that at T2. On the contrary, for 
AT= —4.8 K, the aging at T 2 is no longer independent 
of that at Ti , but gives some contributions to that at T\ . 
As a result, r e ff has a positive finite value. It should be 
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FIG. 9: Difference between e' and e' le f observed by the TC mode for PMMA films with d=514nm and /=100Hz. The aging 
temperature of the first stage Ti is fixed to be 374.7K and T 2 (= Ti + AT) is changed: AT=-4.7 K (d,h), -9.9 K (c,g), -14.5 
K (b,f), and —19.7 K (a,e). Aging times are ri = T2 = 5 hours. In the right figures, r c g is the same as defined in Fig. 8. 



noted here that, although the above result is obtained 
for e", almost similar results are obtained also for e' . 

Next, we show the result obtained for PMMA films 
with d=514 nm, which is much larger than d=26 nm. 
As an example, the observed results of e' for /=100Hz 
are shown in Fig. 9. It is found from Figs.9(e)-(h) that 
for AT—— 19.7 K, r e s is zero, that is, the aging dynam- 
ics at T\ are totally independent of that at T2. How- 
ever, as I AT I decreases, r c g increases monotonically. For 
AT=— 4.7 K, r c ff w T2- This means that the aging at T2 
has almost full contributions to that at T%. On the basis 



of these results, it can be expected that the contributions 
of the aging at Ti to that at Ti change continuously with 
decreasing |AT|. A similar continuous change in r c ff with 
respect to |AT| has also been observed in some examples 
of spin glasses Comparing Fig. 8 with Fig. 9, it 

is found that thinner PMMA films are more sensitive to 
the change in |AT| than thicker PMMA films. 

To measure the contributions of the aging at T% to that 
at Ti , we evaluate the ratio T e ff /T2 as functions of AT un- 
der various conditions. In the case of d — 514nm, we use 
the values of e' for / = 100Hz and 20Hz, and those of e" 
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FIG. 10: Dependence of effective time T e g normalized with 
the aging time of the second stage T2 on the difference of 
two aging temperatures for PMMA thin films with d=26nm, 
514nm and /=20Hz, 100Hz. The effective time r e ff is evalu- 
ated from the data for e and e" . Open and full symbols dis- 
play the results for d=514nm and d=26nm, respectively. The 
curves are obtained by fitting the data to Eq.(8). The solid 
curve is for d = 514nm and the dotted curve is for d — 26nm. 



for / = 100Hz, while in the case of d = 26nm, we use the 
values of e' for / = 20Hz and those of e" for / = 20Hz. 
Figure 10 displays that the effective time r e g normalized 
with T2 is a decreasing function of |AT| and decays to 
zero as |AT| increases. The decaying rate of t s/t2 with 
respect to |AT| is larger for e?=26nm than for d=514 nm. 
Therefore, we can say that the independence of the aging 
dynamics between two different temperatures is accom- 
plished faster in thinner films with increasing |AT|. 

It is interesting to compare the AT dependence of 
T g with that estimated through a thermal activation 
process [2(j. If it is assumed that the aging involves 
a thermally activated jump over a free energy barrier 
At/, the corresponding time r at temperature T char- 
acteristic of the jump motion over the barrier is given 
by t(T) — T<jCxp(AU/kBT), where To is a microscopic 
time scale and ks is the Boltzmann constant. Here, we 
consider the aging at T 2 for a period of r%. In this case, 
we can say that the aging process proceeds by a reduced 
time f(= 72/t(T2)). If the reduced time f of a state is 
the same as that of another state, we can consider that 
both the states are the same. Because the aging process 
at T 2 for a period of T2 corresponds to that at T\ for a pe- 
riod of T ff, we obtain the relation t c r/t(T\) = t 2 /t{T 2 ). 
If the energy barrier AU is independent of T, we obtain 
the relation (t(Ti)/t ) Ti = (t(T 2 )/t ) T2 for the ther- 
mal activation process. Combining the above two equa- 
tions, we obtain the relation for the temperature shift 
from T 2 (=T 1 + AT) to T x : 



T 2 



r(T 2 ) 



AT/Ti 



(8) 



Eq.(8) can reproduce the observed data in Fig. 10. Never- 
theless, from the qualitative comparison of the observed 
results with Eq.(8), it is found that the microscopic time, 
To, associated with the aging is smaller in thin film ge- 
ometry than in the bulk films. 



V. DISCUSSIONS 

A. Physical origin of t w dependence of e' and e" 

As shown in Sec. Ill, we have observed the interesting 
behavior of the aging dynamics. In this section we dis- 
cuss the physical origin of the t w dependence of e' and e" 
observed in the present measurements. The dynamic di- 
electric constant, e*(T), is obtained by dividing the com- 
plex electric capacitance C* (T) at the temperature T by 
the geometrical capacitance Co (To) evaluated at To=273 
K. 

The evaluation of e* in this paper is based on the as- 
sumption that the T and t w dependence of Co can be ne- 
glected compared with that of the intrinsic dielectric con- 
stant. This assumption may not be valid in some cases. 
The T dependence of Co can be neglected as long as the 
value normalized with that at t w =0 is used for a given 
temperature. However, the t w dependence of Co cannot 
be neglected for the following two reasons. Firstly, it is 
well-known that the physical aging of polymers increases 
the density with aging time. This is called densification, 
which leads to the decrease in film thickness with increas- 
ing t w . This contraction should be commonly observed 
independent of film thickness. In addition to this, recent 
measurements done by Miyazaki et al. show that the film 
thickness of polystyrene (PS) thin films supported on Si 
substrate decreases very slowly with time at an annealing 
temperature above T g when the initial thickness is less 
than 20nm [27L |28| . This indicates that there should be 
a very slow relaxation process of PS in thin film geome- 
try. If we take into account the above two processes for 
PMMA thin films, it follows that the geometrical capac- 
itance Co depends on the aging time t w and that the t w 
dependence of Co of ultra thin films can be different from 
that of bulk samples. 

The effect due to the t w dependence of Co is evaluated 
here. Since C = e'eo-f > the relative change in C is given 
by 



AC 



Ae' 



Ad 



(9) 



Because the determination of T e g is accompanied with the 
relatively large uncertainty, it is difficult to judge whether 



where S is assumed not to change with time, and the 
symbol AX(X=C, e', d) is AX = X(t w ) - X(Q). The 
discussions in this paragraph can also be valid for e" and 
C". If the sample goes towards an equilibrium state dur- 
ing the aging process at a temperature, the densification 
proceeds and the thermal fluctuations are suppressed. 
The former leads to the decrease in d, while the latter 
leads to the decrease in e'. In this case, it is reasonable 
to expect that Ad/d and Ae'/e' has the same form with 
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respect to t w . Although the aging strength may be mod- 
ified, we can evaluate the t w dependence of e' through 
that of C . In thin film geometry, however, the t w depen- 
dence of e' obtained thus may be affected by a possible 
existence of the very slow relaxation process of the film 
thickness. 

As shown in Fig. 1(b), we observed that e'(= 
C (T) I 'Cq(To)) increases with t w after decaying up to 
about 11 h during the aging at T a = 375.4K for PMMA 
thin films with the initial thickness of 20nm. This in- 
crease may be due to the existence of the very slow re- 
laxation process of the thickness observed only in thin 
film geometry p7j . The characteristic time of this very 
slow relaxation process should be longer than that of the 
relaxation of the intrinsic e' (e") and the contraction due 
to the physical aging. 

Now let us try to estimate the effect of the existence 
of this very slow relaxation process on the present ex- 
perimental results. Comparing the results in Fig. 1(a) 
and Fig. 1(b), we find that Ae" (t w ) / e" (0) ~ 0.062 and 
Ae'(t w )/e'(0) ~ 0.0048 for t w = llh at T a = 375.4K. 
An analysis of Ad/d in the present measurements shows 
that Ad(t w )/d(0) ~ 0.0009 for t w = llh at T a = 375.4K. 
Therefore, even if the increase in e exists under some 
conditions, the effect of this increase on the dynamics 
of the intrinsic e' is less than 20% when t w is limited 
to t w < llh. It should be noted here that this effect 
can only be observed for e' near T g and for ultrathin 
films (d = 20, 26nm). On the basis of this discussion, we 
believe that the observed dynamics of e' and e" in the 
present measurements can safely be regarded as those of 
the intrinsic dielectric constants for PMMA thin films, 
except for some special conditions mentioned above. 

The discussions here are based on the existence of the 
very slow relaxation mode in PMMA thin films. In or- 
der to check the validity of this assumption, we should 
measure the frequency spectrum of dynamics dielectric 
constant in the lower frequency region. 



B. Comparison with the aging scenario in spin 
glasses 

In the CR and TC modes, we observed the memory and 
rejuvenation effects. There are many common properties 
between polymer glasses and spin glasses. For spin glass 
systems, the memory and rejuvenation effects have been 
discussed by two different models : one is a hierarchical 
model 13, Ijl and the other is a droplet model H H HJ 

The hierarchical model is based on the mean-field 
picture originating from the Sherrington-Kirkpatrick 
model [33| ■ It is assumed that the spin glass (SG) phase 
is characterized with a multi valley structure of the free 
energy surface at a given temperature. At a temperature 
Xi, the system relaxes over the many valleys formed at 
that temperature. When the system is cooled from T% to 
Ti + AT (AT < 0), each valley of the free energy splits 



into new smaller subvalleys. If |AT| is large enough, the 
energy barriers separating the initial valleys are too high 
to go among different initial valleys within the period of 
T2 at T\ + AT. Only relaxations within the initial valleys 
can be activated, and hence the occupation number of 
each initial valleys holds during the aging at T% + AT. 
When the system is heated again from T\ + AT to T\ , the 
small subvalleys merge back into their ancestors. There- 
fore, the memory at the end of the first isothermal aging 
process at T\ can be memorized and recalled at t\ + t%. 
The rejuvenation effects can also be explained in a similar 
way. 

As for the droplet model, the growth of SG ordered 
domain is taken into account. In a spin glass system, the 
growth rate of the domain is very slowed down because 
of the heterogeneous distribution of the frustration. Fur- 
thermore, the memory and rejuvenation effects can be 
explained within the droplet theory using the tempera- 
ture chaos concept, in which the spin configuration or 
the energy landscape can be changed globally even for 
any infinitesimal temperature change on a length scale 
larger than the overlap length. 

Recent studies based on the droplet model produced 
more powerful scenarios such as the droplets-in-domain 
scenario Q and the ghost domain scenario j3j||3(|. Fur- 
thermore, disordered systems showing the memory effects 
have been classified into 2 categories as follows: 

1. the system showing strong rejuvenation effects and 
memory effects: The strong rejuvenation is due to 
the temperature chaos effect. The relaxation asso- 
ciated with the strong rejuvenation after the tem- 
perature shift from T\ to T2 is the same as that 
observed after the direct cooling to the aging tem- 
perature T2. 

2. the system showing no strong rejuvenation, but a 
transient relaxation and memory effects: In this 
case, the strength of the transient relaxation after 
the temperature shift by AT becomes stronger with 
increasing |AT|. The aging proceeds cumulatively 
even if the temperature is shifted. 

In some systems, both properties can be observed de- 
pending on the temperature shift AT. For a sufficiently 
small AT the aging proceeds cumulatively. When AT 
becomes large, the aging is no longer cumulative but the 
noncumulative and strong rejuvenation is observed just 
after the temperature shift. 

Polymeric systems have been so far classified into 
the type 2). Although the memory effects have been 
observed, it has been reported in the literature that 
no strong rejuvenation occurs after the temperature 
shift [35|. However, the present measurements show that, 
after the negative temperature shift from T\ to T2 by 
AT ps — 20K, the relaxation of e" due to the aging is reini- 
tialized, i.e., goes back to the value on the reference curve 
measured during the cooling process. The values of the 
e", then, begin to decay at T2 along the relaxation curve 
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measured after the cooling directly to T2, as shown in 
Fig. 7. From this results, we believe that there is a strong 
rejuvenation in the case of PMMA for AT w -20K. Fur- 
thermore, as AT decreases from — 4.8K to — 19.6K in 
Fig. 8, the contributions of the aging at X2 to that at T\ 
decrease, and the effective time t c s decreases to with 
increasing |AT| in the negative temperature cycle. This 
experimental results imply that there is a crossover from 
the cumulative dynamics to the noncumulative ones as 
the value of |AT| increases. In order to check whether 
there is really the crossover from cumulative to noncumu- 
lative, i.e., whether there is the temperature chaos effect 
in polymeric systems, the 'twin-temperature shift' exper- 
iments should be performed, which have been developed 
in spin glasses [371 ]. 

Comparing polymeric systems with spin glasses, the 
characteristic time of the microscopic flip motion of 
dipoles or spins can be expected to be much larger in 
polymeric systems than in spin glasses. Therefore, poly- 
meric systems may be a model system for investigating 
the aging dynamics in the shorter time regime, since an 
experimental time window corresponds to shorter time 
scales in polymeric systems compared with atomic spin 
glasses. Superspin glasses (strongly interacting nanopar- 
ticle systems) [23 can also be regarded as a model system 
for the same purpose, because they have the longer mi- 
croscopic flip time of a superspin than atomic spin glasses 
have. The polymeric systems may have a still longer mi- 
croscopic flip time, and hence we expect that the poly- 
meric systems can be better model systems. 

C. Thickness dependence of aging dynamics 

As for the thickness dependence of the aging phenom- 
ena, we could obtain several interesting results: 1) the re- 
laxation strength due to the aging decreases with decreas- 
ing film thickness. 2) the width of the dip in e' — e' lci dur- 
ing the cooling process in the CR mode becomes smaller 
with decreasing film thickness. 3) the effective time T e g 
in the TC mode decays to zero more quickly in a thin 
film geometry with increasing |AT| in the negative AT 
region than in bulk systems. 

As mentioned in Sec. IV, some of the above results re- 
lating to the thickness dependence may be explained by 
assuming that the fundamental microscopic flip time To 
of dipoles becomes smaller with decreasing film thick- 
ness. As shown in Fig. 5, the width of the dip observed 
during the cooling process in the CR mode is ~ 20K 
for d = 514nm and ~ 10K for d — 26nm in the case of 
PMMA films. In (atomic) spin glasses, where To is much 
smaller than in polymer glas ses, the width of the dip of 
e' in the CR mode is ~ 3K |2q|. which is even smaller 
than the value for the thinnest film thickness. The width 
of the dip observed in the CR mode does not depend 
on the cooling rate for the range of cooling rate from 
0.08K/min to 0.33K/min, although the depth of the dip 
strongly depends on the cooling rate |Tcj . 

VI. SUMMARY 



dielectric measurements. The results obtained in the 
present measurements are as follows. 

1. The relaxation of dielectric constants is observed 
during the aging process at an aging temperature 
T a . The relaxation strength decreases with decreas- 
ing T a and with decreasing film thickness. In the 
films of d = 20nm and 26nm, the increase in e' has 
been observed at T a = 375K, which may be due to 
the existence of the ultra slow relaxation process 
of the film thickness observed recently through the 
X-ray reflectivity measurements. 

2. In the CR mode, the aging at T a by the way of 
cooling to room temperature can be memorized as 
the dip in e' — ej. cf curve around T a . During the 
subsequent heating process the memory can be re- 
called. The width of the dip in e' — e' Te{ decreases 
with decreasing film thickness. 

3. In the TC mode with negative AT, a strong re- 
juvenation effect of the difference between e" and 
e" e{ has been observed just after the negative shift 
from Ti to T 2 by AT » -20K. At the same time, 
a complete memory effect could also be observed 
when the temperature goes back from T2 to T\ . 

4. As |AT| decreases in the negative AT region in the 
TC mode, the contributions of the aging at T2 to 
that at Ti become larger i.e., the effective time T e ff, 
which is a measure of the contribution, increases 
from to T2- A continuous change in r e ff with AT 
could be observed, which reminds us of a crossover 
between the noncumulative aging and the cumula- 
tive one. The dependence of r e ff on AT depends 
on the film thickness. 

In this paper, we focused on the aging dynamics ob- 
served in dielectric constants. However, similar aging 
phenomena can also be observed in heat capacit y o f 
PMMA observed by temperature modulated DSC |38|. 
It can be expected that the aging phenomena observed 
in this paper are common ones for any susceptibility in- 
cluding dielectric constant, magnetic susceptibility, and 
heat capacity. 
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